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Recent reports on intense-field pump-probe experiments for liigh harmonic generation from coher- 
ently rotating linear molecules, have revealed remarkable characteristic effects of the simultaneous 
variation of the polarization geometry and the time delay on the high harmonic signals. We analyze 
the effects and give a unified theoretical account of the experimental observations. 
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The phenomenon of high harmonic generation (HHG) 
from atoms or molecules in intense laser fields can be 
thought of as a "fusion" of n laser photons, each of en- 
ergy huj, into a single harmonic photon of an enhanced 
energy Ml = nhio. This might seem surprising at first 
since the photons do not interact with each other and 
therefore can not "fuse" on their own. However, a bound 
electron interacting with a laser pulse can absorb n pho- 
tons from the laser field, go into highly excited virtual 
states and can return to the same bound state by re- 
leasing precisely the excitation energy (nhui) as a single 
harmonic photon. Note that at the end of the coherent 
process the electron does not change its state - it merely 
acts as a "catalyst" of the process. The phenomenon 
is currently being vigorously investigated, both experi- 
mentally and theoretically, specially in connection with 
dynamic alignments of molecules (e.g. [ij). 

Recently a number of remarkable pump-probe exper- 
iments for high harmonic generation with intense fem- 
tosecond laser pulses from coherently rotating linear 
molecules (e.g. N2, O2, CO2, HC=CH) have been re- 
ported in this journal and elsewhere (e.g. 0, [E S S [3l ) ■ 
These experiments measure the HHG signals as a func- 
tion of the delay-time, t^, between a pump pulse that 
sets the molecule in coherent rotation, and a probe pulse 
that generates the high harmonic signal from the rotating 
molecule. The changes in the dynamic signals are then 
recorded by varying the angle, a, between the polariza- 
tions of the two pulses. Fig. [l]shows a schematic diagram 
of the various vectors involved in the pump-probe experi- 
ments. The geometric angle a is the operational angle in 
the laboratory, although at times it is erroneously identi- 
fied with the angle 6 (or 6'); the latter is a quantum vari- 
able, not measured in these experiments. Here we derive 
an explicit theoretical expression for the HHG signal as 
a simultaneous function of the geometric angle a and the 
delay-time td and analyze the experimental observations. 
The results provide a unified theoretical account of the 
observed effects. 

Let the total Hamiltonian of the molecular system in- 
teracting with a pump pulse Li at a time t, and a probe 
pulse L2 applied after a delay-time td, be written, within 



Figure f : A schematic diagram defining: molecular axis, R, 
electron position r, pump polarization ei, probe polarization 
62; z and z' axes lie on common z-z'-x plane; fields propagate 
along y-axis. 



the Born-Oppenheimer approximation as (e.g. [1, [^): 

Htot it) = ff^"' + Vn^L, (t) + iff^ + Ve-L, (t - td) (1) 

where the subscripts N and e stand for the nuclear and 
the electronic subsystems, respectively. An intense fem- 
tosecond pump-pulse is assumed to interact with the 
molecular polarizability, via VN-Li(t), and sets it into 
coherent free rotation. 

The coherent rotational motion 10] is described by the 
nuclear wavepacket states created by the pump pulse: 
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Each wavepacket state ^ evolves one-to-one from an 
initially occupied ensemble of eigen states, \JoMq), 
populated with a Boltzmann distribution p{Jo) = 
^ where Zp is the partition function. Thus, 
after the pump pulse, the initial state of the molecule is 
characterized by the ensemble of product states, |xi(0)i 
with i = {e, Jo, -^0}, composed of the ground electronic 
{t) \ and the coherent wavepackets \^j„Moit)}' 
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Generalizing the well-known strong-field KFR (Keldysh- 
Faisal-Reiss) approximation (e.g. [lH) to the present 
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Figure 2: Quantum amplitude for coherent emission of a 
liigh harmonic photon (frequency fl) is the sum of two di- 
agrams, (a) direct, and (b) time- reversed; probe-interaction 
(line-x), photon emission (arrow); intermediate propagators, 
Go; Volkov wave- vector k; i = Eq.Q. 



molecular case, we write the wavefunction of the system, 
evolving from each of the ensemble of the initial states 
(13 as: 

= IX^W) + / dt'Go{t,t')Ve-L2{t' -td) \x^{t')) 



(4) 

where, the Green's function Go{t,t') of the system is 
given by 
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where. 



are ionic orbitals and \<j}p{t)) are Volkov 

states (e.g. [llj). 

The quantum transition amplitude for the coherent 
emission of a harmonic photon of energy Ml = nhuj, from 
an initial state ([3]) evolving into ([4]) and recombining back 
into the same state (jS]), is given by the sum of a 'direct' 
and a 'time reversed' diagram for the photon emission 
process (cf. Fig. [2]). Writing out the amplitude ana- 
lytically using Eqs. ([2]) to dS]), assuming the "adiabatic 
nuclei" condition, Max{AEjj') ^ E,,, carrying out the 
lengthy algebra, and modulo-squaring the result, we ob- 
tain the coherent HHG emission probability for each ini- 
tial state (|31). Taking the statistical average of the in- 
dependent probabilities for the ensemble of initial states 
^ , we obtain the scaled HHG signal "per molecule" , as 
an explicit function of a and td- 
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with, L = {\l - l'\, (l + I')), M = {-L,L); the parame- 
ters a^^!'l, {I , I' , L; m) are given by rather lengthy but ex- 
plicit expressions [Qj] that depend on the partial angular 
momenta l{l') of the active electron and their conserved 
projection, m, on the molecular axis, on the matrix el- 
ements of the absorption and recombination transition- 
dipoles, and on the usual vector addition coefficients. 

Specializing Eq. ([7]) to the case of N2 (molecular or- 
bital symmetry (Tg, m = 0; dominant 1(1') — 0,2,4), we 
get, in an ordinary trigonometric representation, an an- 
alytic expression of the dynamic HHG signal for N2: 



{td,< 



P1+P2((C0S2 e')){td) 
+P3((C0S2 e'){cOS^ 0')){td) 

+p^{{cQs'e')){td) + --- 

+pio{{cos''9'){cos^d')){td) 



(8) 



where cos 6' — cos a cos 9 + sin a sin 9 cos (p. Similarly, for 
O2, [t^q symmetry, m — 1, and dominant — 2,4) we 
get. 
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<J(sin2 9' cos2 9')^) (td) 
q2 (^(sin^ 9' cos^ d') 
(sin2 e'cos'^ 9')'^{td) 

{^(sin^ 9' cos^ {td) 



(9) 



The coefficients p's and q's are determined by simple com- 
binations of the parameters a^?], (Z, I', L; m) [9|]. We note 
that for a =0,Eqs. ([8]) and ([9]) reduce correctly to the 
special limits [8|. 

In Fig. [3] we show the results of computations using 
Eq. ([8]) for the dynamic signals from N2 as a function 
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Figure 3: Calculated 19th harmonic dynamic signal for N2 
for various pump-probe polarization angles, i.e. a = 0°, a = 
45°, and a = 90°; pump intensity / = 0.8 x 10^''W/cm^ 
probe intensity, / = 1.7 x 10^'* W/cm^, duration 40 fs, and 
wavelength 800 nm; Boltzmann temperature 200 K. 
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of the delay time td, at three different relative polar- 
ization angles, a =0°, 45°, and 90°. The results show 
the full revival with a period T^ev = jSc ^ 8.4 ps, and 
the fractional i and j revivals, for all three a values. 
The ((cos^6')) {td) term is known to govern the ^ and 
J revivals and the associated Raman allowed spectral 
lines (e.g.P, @]). Remarkably, the signals for a —0° and 
a = 90° are found to be in opposite phase, while that for 
a = 0° and a — 45° are in the same phase. Exactly the 
same phase relation between the a-dependence of the id- 
signal from N2 has been observed in recent experiments 
(e.g. iiH). 

To analyze their origin, we consider the leading term 
of Eq. ([5]) for N2, more explicitly. (Below, we omit 
the argument {td) for the sake of brevity.) Noting that 
(cos^ 6') — \ siv? a + (cos^ ~ \ ^in^ a) (cos^ 0) we get: 

S{td;a) ~ ^pi +]32^ sin^ +P2 ^cos^ a - i sin^ 

x((cos2e')) + --- (10) 

Therefore, for the parallel polarizations we have, 
S {td] 0°) ~ pi + P2 ((cos^ 0)) and for the perpendicular 
polarization, S{td;90°) w {pi + ■^P2) - 5P2 ((cos^ 
Clearly due to the opposite sign of the ((cos^ 9)) term, 
they vary in opposite phase to each other from their 
respective bases. In contrast, the signal at a = 45°, 
S'(trf,45°) w {pi + JP2) + jP2 ((cos^e)), has the same 
sign of the ((cos^ 0)) term as for a —0, which makes them 
to vary in phase. These behaviors are what can be seen 
in the full calculations in Fig. 3, and they also agree with 
the recent experimental observations (e.g. 0, H, 0, Bl). 
The simple formula ([10]) predicts further that the ex- 
trema of the signal should occur for sin a cos a = 0, with 
a maximum at a =0° and a minimum at a = 90°. This 
is also what has been seen experimentally 0, 0| • Finally, 

pn)) predicts a "magic angle" ac = arcsiny'^ « 54.7°, 

given by the condition (cos^ Q^c — ^ sin^ ac) = at which 
the HHG signals become essentially independent of the 
delay td between the pulses. Exactly such a "magic" 
crossing angle for N2 signals has been observed experi- 
mentally . We may point out that this geometry can be 
used in femtosecond pulse-probe experiments to generate 
an essentially steady HHG signal from freely rotating N2 . 

In Fig. [4] we present the results of full calculations for 
O2, using Eq. ([9]), for the three geometries, a = 0°, 45°, 
and 90°. The signals are seen to be characterized by a full 
revival at Trev = ~ 11-6 ps and also by the fractional 
i and J revivals, like in N2, as well as an additional ^- 
revival, for all the three geometries; the same characteris- 
tics have been observed experimentally (e.g. 0, 0|)- The 
existence of the ^-revival is due mainly to the presence of 
higher powers and moments than ((cos^ 9)), that couple 
the Raman-forbidden (AJ = ±4) and the "anomalous" 
transitions (|AJ| > 4) between the rotational states Q. 
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Figure 4: Calculated 19th HHG spectrum of O2 for various 
pump-probe polarizations angle, i.e. a — 0°, a = 45°, and 
a — 90°; pump intensity / = 0.5 x 10^''W/cm^, probe inten- 
sity, J = 1.2 X 10^'' W/cm^, duration 40 fs, and wavelength 
800 nm; Boltzmann temperature 200 K. 

We may express the contribution from the first term of 
Eq. ^ more explicitly as: 

S{td;a) « 7-7 ( {(3 — 30cos^ a + 35 cos^ a) 
64 \ 

X (sin^ 9 cos^ 9) 

— (1 — 6 cos^ a + 5 cos^ a) (cos^ 9) 

+ (4sin2a - 3sin'*a)}^^ + • • •. (11) 

For a — 0°, this gives, S{td;a) ~ 
gi {^(sin^Scos^ei)^^ and, for a = 90°, 5 (t^; a) « 

g^(|)2 ^{(sin^^cos^e*) - i(cos2 6i) + ij^^ A compar- 
ison of the above expressions suggests that for a = 0° 
and 90°, the signals at the full, ^ and j revivals would 
be in opposite phase, and that at the ^ revival would be 
in phase. A direct comparison of the calculations using 
the above abbreviated formulas with the full calculation 
in Fig. 4 and the experimental observations in O2 (e.g. 
11,01), fully confirm the above expectations. 

In Fig. [5] we show the calculated results of the dynam- 
ical signals for N2 (upper panel) and O2 (lower panel), 
as a continuous function of a, between 0° to 90°, at 
three different delay-times td near the ^ revival period. 
In the upper panel for N2, a remarkable coincidence of 
the three signals is seen to occur at the "magic angle" 

Q!c — arcsiny^ « 54.7° as predicted above. More- 
over, the signal at the "top" -alignment time td — 4.05 ps 
(solid curve) is seen to lie above the signal at the "anti- 
top" alignment time td = 4.3 ps (dash-dot curve), for all 
a < ac, and they invert their relative strengths for all 
a > ac- This is again in agreement with the recent ob- 
servations (e.g. [3, 0]). The corresponding signals for 
O2 (lower panel) does not show a single crossing point. 
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Figure 5: Variation of dynamic HHG signal as a function of 
pump-probe angle a, near the first half-revival, for and 
O2. The pulse parameters are the same as in Fig. ([Sj for 
and Fig. © for O2 ■ 

rather they cross at three different points a, b, and c in 
the neighborhood of the magic angle ac ~ 54.7°. Such a 
crossover-neighborhood around the magic angle « 54.7° 
for O2 is recently confirmed experimentally [5[. The ab- 
sence of a single crossing point for O2 is due mainly to 
the non- negligible contribution of the moment ((cos'' 0)) 
to the O2 signal (cf. Eq. (dJ)). 

Before concluding, we may make a few qualitative re- 
marks on the a dependence of the dynamic signals for 
the more complex triatomic molecule CO2 Q , and the or- 
ganic molecule acetylene (tt symmetry) , that are mesured 
recently [1,0|- The structure of the operator ([7]) sho-ws, 
even without a detailed calculation, that the CO2 and 
acetylene (HC=CH) , due to their linear structure, would 
show a similar crossover at or near « 54.7°. A direct 
perusal of the experimental data 0, Q confirms this gen- 
eral expectation from the present theory - both CO2 and 
acetylene exhibit the crossover effect, and indeed near the 

"magic angle" ac — arcsin Ki 54.7°. 

To summarize: The simultaneous dependence of the 
dynamic HHG signals from coherently rotating linear 
molecules, on the relative polarization angle, a, and the 
time delay, td, between a pump and a probe pulse, is 
investigated theoretically. A general formula for the dy- 
namic signals for linear molecules is derived (Eqs. ((6]) and 
([7|)). It is used to analyze the recently observed charac- 



teristics of the HHG signals from N2 and O2. Among 
other things, a "magic angle" ac = arcsin ~ 54.7° 
for the crossing of the dynamic signals for N2, and a 
crossover-neighborhood around the "magic angle" , for 
O2, are predicted by the theory and confirmed by the 
available experimental data. The presence of analogous 
crossovers for the more complex linear molecules, CO2, 
HC=CH (acetylene), are also suggested by the present 
theory, and are corroborated by the recent observations. 

We thank Prof. K. Miyazaki for the private communi- 
cations and fruitful discussions. This work was partially 
supported by NSF through a grant for ITAMP at Har- 
vard University and Smithsonian Astrophysical Observa- 
tory. 



[1] F. Rosca-Pruna and M.J.J. Vrakking, Phys. Rev. Lett. 
87, 153902 (2001); R. Velotta et al, Phys. Rev. Lett. 87, 
183901 (2001); N. Hay et al., Phys. Rev. A 65, 053805 
(2002); N. Hay et al., J. Mod. Opt. 50, 561 (2003); P.W. 
Dooley et al., Phys. Rev. A 68, 023406 (2003); M. Lein 
et al., Phys. Rev. A 67, 023819) (2003); J. Itatani et al., 
Nature 432, 867 (2004); R. de Nalda et al., Phys. Rev. 
A 69, 031804 (2004); M. Kaku et al., Japan j. Appl. 
Phys. 43, L591 (2004). Zeidler et al. in Ultrafast Optics, 
IV, ed. F. Krauz et al.,( Springer, New York, 2004); M. 
Lein et al., J. Mod. Opt. 52, 465 (2005); C. Vozzi et al., 
Phys. Rev. Lett, 95, 153902 (2005); X.X. Zhou et al., 
Phys. Rev. A 72, 033412 (2005); C.B. Madsen and L.B. 
Madsen, Phys. Rev. A 74,023403(2006). C. Vozzi et al., 
J.Phys. B. 39, S457 (2006). 

[2] J. Itatani et al., Phys. Rev. Lett. 94, 123902 (2005). 

[3] K. Miyazaki et al., Phys. Rev. Lett. 95, 243903 (2005). 

[4] T. Kauai et al.. Nature 435, 470 (2005). 

[5] K. Miyazaki et al., IEEE-CLEO-PR-2007 Conf. Rep. 
(submitted, 2007). 

[6] N. Kajumba et al.. Central Laser Facility Ann. Rep., 
2005/2006, p. 80, Rutherford Appleton Lab. U.K. (2006). 

[7] R. Torres et al., Phys. Rev. Lett., 98, 203007 (2007). 

[8] F.H.M. Faisal, A. Abdurrouf, G. Miyaji, and K. 
Miyazaki, Phys. Rev. Lett., 98, 143001 (2007). 

[9] A. Abdurrouf and F.H.M. Faisal, Phys. Rev. A (to be 
submitted, 2007). 
[10] H. Stapelfeldt and T. Seideman, Rev. Mod. Phys. 75, 
543 (2003). 

[11] A. Becker and F.H.M. Faisal, J. Phys. B 38, Rl (2005). 

[12] G. Herzberg, Molecular Spectra and Molecular Structure, 
II (Van Nostard Reinhold, New York, (1950). 



